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& ■ Abstract 

OJQ. 

A central question of dynamics, largely open in the quantum case, is to what extent it erases a 

system's memory of its initial properties. Here we present a simple statistically solvable quantum 



model describing this memory loss across an integrability-chaos transition under a perturbation 
obeying no selection rules. From the perspective of quantum localization-delocalization on the 



lattice of quantum numbers, we are dealing with a situation where every lattice site is coupled to 
every other site with the same strength, on average. The model also rigorously justifies a similar 



set of relationships recently proposed in the context of two short-range-interacting ultracold atoms 



J> , in a harmonic waveguide. Application of our model to an ensemble of uncorrelated impurities on 

<N 

C ' a rectangular lattice gives good agreement with ab initio numerics. 
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There are two basic prototypes of non-dissipative dynamics, corresponding to whether 
time propagation retains a very strong or a very weak memory of the initial conditions: 
integrable dynamics, where time evolution conserves as many independent quantities as there 
are degrees of freedom; and completely chaotic dynamics, where all the traces of the initial 
state of the system — except for its energy and perhaps a very few other quantities — are erased 
exponentially fast. As an integrable system is perturbed away from integrability, its many 
conserved quantities stop being conserved, and a fundamental problem is to understand the 
manner and the mechanism by which that happens. 

In the case of classical systems, this problem was highlighted by the FPU "paradox" and 
was central for the development of several subfields of classical mechanics, including the 
theories of solitons and of dynamical chaos [1] . The most celebrated result is the KAM the- 
orem, which explains the gradual disappearance of quasiperiodic orbits as the perturbation 
increases. 

In quantum mechanics, there do exist several relevant experimental [2|-l4] and numerical 

studies, all in the context of studies of relaxation, as well as attempts to formulate a 
quantum KAM theorem 6|; nevertheless, the problem remains largely open. 

For isolated quantum- chaotic systems, the key property governing the memory of the 
initial state under time propagation is eigenstate thermalization, which is defined as follows: 
consider the values (a|A|a), which are the quantum expectation values of an observable A 
with respect to the eigenstates |a), with the corresponding eigenenergies E a , of the system 
hamiltonian H. The system has the eigenstate thermalization property with respect to the 
observable A if, in any limit in which the density of states increases without bound, and in 
any microcanonical energy window, all the values (a|A|a) in the window become equal to 
each other (here (a|A|a) is said to be in a microcanonical energy window if the eigenenergy 
E a , corresponding to the state |a), lies within some specified narrow energy interval). This 
fact was first discovered in the context of studying the transition from quantum to classical 
behavior, in the initiating paper by Shnirelman [7| and in the body of work in mathematics 
and mathematical physics that immediately followed (see ref. [8J and the references therein). 
Eigenstate thermalization is now believed to be the main explanation for why large quantum 



systems behave thermodynamically [9l4l l| . In the present context, however, of greatest 
interest is the residual variation among the (a|A|a) values within an energy window — a 
consequence of the fact that the density of states is not actually infinite. This variation 



has been linked to a physically relevant autocorrelation function 12[ and studied in many 



examples 



|8|,|l3rll5 



15] ; however, the most relevant point for the present purposes is that the 
residual variation provides an upper limit on how much difference there can be between the 
following two quantities: the first is the infinite time average of the quantum expectation 
value of an observable A with respect to the time-evolving state \ip(t)), 
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where A = ^ J dt (ijj(t)\A\ijj(t)), iV'init.) is the initial state, and \ip(t)) = ^2 a c a (t) \a) with 
c a (t) = (c^lV'imt.) e~ ttEa ^ h ] the second quantity is the microcanonical average of (a|A|a). The 
latter is simply the average of the (a|A|a) values over all states |a) in a narrow energy 
window centered at E = (V'mit.l-^lV'imt.), the mean energy of the initial state. Note that the 
occupation numbers |(a|^i n i t .)| are non- negligible generally only for eigenstates \a) from a 
narrow energy window centered at the mean energy E (see the Supplementary Discussion 



for ref . 



11|). Now, if there is some nonzero residual variation among the values of (a|y4|a) 



within a microcanonical window, then one can pick an initial state whose occupation numbers 
| (ct l^init.) | tend to be, say, larger for states \a) whose values (a|A|a) are larger. Then the 
infinite time average A will be larger than the microcanonical average of the (a|A|a) values. 



In contrast, as pointed out in refs. |9l-ll 1| . if all the values (a|A|a) are the same, then A 
and the microcanonical average of the (a|/l|a) values are the same. 

In the other extreme, in the case of integrable systems, time evolution preserves the 
values of all quantum numbers [2|, although it has been shown that, at least in the case of 



large systems, the infinite time average is able to erase a 



It 



re information about the relative 



phases of the eigenstates composing the initial state 16Ml8| . Eigenstate thermalization fails 
for integrable systems. 

A real micro- or mesoscopic physical system is usually neither completely integrable nor 
completely chaotic. On the other hand, it may often be described as an integrable system 
perturbed by an integrability-breaking potential. It is therefore of central importance to 
understand how the dynamics of the system changes as the strength of the perturbing 
potential is increased from zero to the point where all resemblance to the original integrable 
system is lost. 

In the present work, we show that the gradual disappearance of the integrals of mo- 
tion under perturbation can, for a conceptually important class of quantum systems, be 



completely characterized by a simple and universal relation. Our approach is to develop 
an exactly solvable statistical model that captures the essence of the physics, and then 
test its predictions — in the full range of parameters from the integrable regime through the 
well-developed chaos — against the exact time dynamics of a Deformed Random Gaussian 
Model[l9j, a Seba-type billiard 20|, and a two-dimensional Anderson model 2lH23j]. All 
three are examples of "rough billiards," by which we mean that in their interiors they have 
scatterers that are completely decorrelated, and, as a result, the perturbation has large ma- 
trix elements between any two momentum states. Correspondingly, in our statistical model, 
we will be considering perturbing potentials that obey no selection rules. From the point of 
view of "localization-delocalization in the space of quantum numbers" our model leads to a 
lattice with extremely long-range hoppings. On the other end of the spectrum one finds the 



kicked rotor 24| and the Bunimovich stadium [25| . with a Coulomb scatterer in a box as an 



intermediate case [26j. In the thermodynamic limit, our model always leads to delocalized 
eigenstates: the phenomena described in this article are thus mesoscopic finite-size effects. 
We show that in our model, the memory of the initial conditions is describable in terms of 
a simple analytic expression that is valid through the full range of perturbation strengths, 
from the integrable to the completely quantum-chaotic extremes, and that features a great 
deal of universality. This type of universality in the memory of the initial conditions in 
quantum systems was first noticed in the context of two short-range-interacting ultracold 
atoms in a harmonic waveguide [27J . In the present paper, we give another physically rel- 
evant example of a perturbation with no selection rules: a lattice billiard with uncorrelated 
impurities. Here the absence of long-range correlation allows the perturbation to couple any 
momentum state with any other, with, on average, the same strength. 

RESULTS 

Formulation of the problem and intuitive analysis We consider the hamiltonian of an 
integrable quantum system with an integrability-breaking perturbation V: H = H Q (n) + V. 
Here the integrable part of the hamiltonian is a function H of all the members of the 
complete set of the integrals of motion, where the latter are labeled asn= (hi, n 2 , . . . , fid)', d 
is the number of the degrees of freedom. The eigenstates and eigenvalues of Ho(n) are \n) and 
Eft, respectively; the corresponding quantities for H are \a) and E a , as before. We assume 



that both the Eft and the E a spectra are free of degeneracies. The non-degenerate spectrum 
is a generic property of an integrable hamiltonian unless it has mutually non-commuting 



integrals of motion or commensurate frequencies [28l |29j. However, because of the absence 
of level repulsion, the energy levels are often near-degenerate: they are distributed as if by 
a Poisson process, resulting in the exponential distribution of level spacings [30| • 

Now we imagine that the system is prepared in some initial nonequilibrium state l^init.) = 
\i/j(t = 0)), and then allowed to evolve according to the hamiltonian H. Following a transient 
period, the expectation value of a generic observable A will, overwhelming majority of the 
time, fluctuate around the infinite-time average given in equation ([T]). Figure [T] gives an 
example of such a process, in the regime intermediate between fully integrable and fully 
chaotic, showing a partial retention of the memory of the initial state. 

Since we are concentrating on the disappearance of the integrals of motion, we take our 
observable of interest to be diagonal in the eigenbasis of H (n), the integrable part of the 
hamiltonian: (n\A\n') = AfiSftft'- Let us choose the initial state to be an eigenstate, |ra init .), 
of Ho(n). In the intuitive discussion that follows, it will help to imagine that we chose a state 
with a "highly improbable" (meaning: very different from the microcanonical average) value 
of the observable of interest A. According to equation ([1]), in the course of time evolution, 
the initial state |?? iriit .) gets transformed into a superposition, with essentially random phase 
relationships, of the eigenstates \a) of the perturbed hamiltonian, H, such that each |a) in 
the superposition has a substantial overlap with |rz ini t.) , i.e. a large corresponding occupation 
number. 

Now, since the eigenstates \n) of H (n) (let us call these the "unperturbed eigenstates") 
form a basis, we may expand the eigenstates |a) of the full, perturbed hamiltonian H 
(the "perturbed eigenstates" ) as linear combinations of the unperturbed eigenstates. The 
stronger the perturbing potential V, the more different the perturbed eigenstates from the 
unperturbed ones, and the greater the number of unperturbed eigenstates that must be 
linearly combined (with appreciable weights) to build any perturbed eigenstate. A stan- 
dard measure of the latter number is the so-called number of principal components, to be 
defined momentarily. Among the perturbed eigenstates \a) whose overlap with |n; n it.) is 
non- negligible, a typical one will be such that when it is expanded in terms of the unper- 
turbed eigenstates, the weight of the state |fzi n i t .) is of the order of rja = Y^n I {n\ a ) I > which 
is the so-called Inverse Participation Ratio (IPR) 15]. The number of principal components 




FIG. 1. An example of a partially retained memory of the initial state. Shown are the x- 
and y-components of the hopping energy (respectively: E x in solid red, with initially positive value, 
and E y in dotted blue, with initially negative value), as a function of time, of a two-dimensional 
33 x 33-site noninteracting Anderson model in the presence of an Aharonov-Bohm (A-B) flux (see 
the first subsection of the Methods section). A fully chaotic systems would obey the equiparti- 
tion theorem, predicting that the infinite time averages of the two components should be equal, 



(E x — Ey 



0. However, for intermediate strengths of the integrability-breaking perturbation, a 



residual deviation from equipartition — strongly correlated with (E x — E y ) , the initial deviation — 
remains. The two straight horizontal lines on the right-hand side of the plot are the exact values 
of E x and E y , computed from equation (pQ) . In terms of the Anderson model parameters in- 
troduced in the text following the paragraph containing equation (|3j), this system has e = 1.0, 
corresponding to the Anderson's disorder parameter of W/J = .87. The initial state was repre- 
sented by an eigenstate of the impurity-free lattice, iV'init.) = \n x = — 9, n y = +3), whose energy was 
E = —1.46 J, where J is the hopping constant (see the first subsection of the Methods section). At 
zero A-B flux, the ground state is \n x = 0, n y = 0). The A-B fluxes (j) x and <p y through a complete 
loop along the x- and the y-directions, respectively, were chosen to be far from approximately 
commensurate; we had <p x = (ip/4)(fto and <\> y = (e/lO)0o, where (fio is the elementary quantum 
of magnetic flux, ip = 1.618. . . is the golden ratio, and e = 2.718. . . is the base of the natural 
logarithm. The apparent correlation between the two energies is due to the conservation of the 
total energy. 



mentioned above is then given by (N t 



pc ; 



• ' ~ 1/rja \ in this case, it is the number of principal 



unperturbed components. 

As we said, we will assume that the matrix elements Vn^ of the integrability-breaking 
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perturbation V obey no selection rules. We will also assume that the equi-energy surface 
Ho(n) is "sufficiently irrational," which is a concept that needs a bit of explaining: as is 
well known, energy eigenstates in an integrable system are completely characterized by their 
quantum numbers (which we can take to be integers) corresponding to all the conserved 
quantities. In a typical integrable system, for sufficiently high excited states, if one of the 
quantum numbers is changed by 1 while all the others are held fixed, the change in the 
energy of the eigenstate will, in general, be large compared to the energy level spacing; that 
is, there will be many eigenstates whose energies lie in between the energy of the original 
state and the energy of the state that differs from the original by 1 in a single quantum 
number. Eigenstates that are neighbors in energy generally have very different values of 
all quantum numbers. In fact, if one orders the eigenstates by increasing energy and looks 
at how their quantum numbers change from one eigenstate to the next in the sequence, 
these changes may appear quite random; if they do, we say that the equi-energy surface is 
"sufficiently irrational." More precisely: let {{nk)v} be the sequence of the values of the kth 
quantum number as one is going from one eigenstate of H (n) to the next in the order of 
increasing energy, labeled by v. The equi-energy surface is said to be sufficiently irrational 
if, for every k, the sequence {{fik)v} passes any simple statistical test for randomness. This 
property seems to be quite generic for integrable systems with incommensurable frequencies; 
see Supplementary Figure S2 for an example. 

These two assumptions (the absence of selection rules in the perturbing potential V and 
the sufficient irrationality of the equi-energy surface of H (n)) imply that the presence of 
the state |ninit.) as one of the components of the state |a) does not "bias the selection" of the 
other states \n) that enter the expansion of \a). This means that the expansion of \a) ap- 
pears as if the states (other than |n; n i t .)) that enter into it were indiscriminately chosen from 
a microcanonical shell around |ni n ; t .). This would make — were it not for the "systematic" 
presence of the initial state in the expansion of |a) — the value of A, on average, equal to the 
microcanonical average. The state |ri init .) becomes the only component of the state \a) that 
"remembers" the initial value of A. Since on average (A^c) ~ l/i] unperturbed states 
enter the expansion of a perturbed state — and out of these unperturbed states that enter 
the expansion, usually exactly one has the initial value for the observable A — we obtain the 
estimate that the infinite time average of A is the weighted average of the initial value and 



the microcanonical value of A, where the weight ratio is 1 : ( (A^, c )^™ — 1 ) in favor of the 
microcanonical value: A ~ r)a A +(1 — r]a)A MC , where A = {ni n it.\ A\n init ) is the quan- 
tum expectation value of A with respect to the initial state, and Ayic is the microcanonical 
average, i.e. the average over all the values (n\A\n) such that the unperturbed eigenenergies 
of \n) lie in a narrow window centered at the mean energy of the system. 

A statistically solvable model. Guided by this intuitive reasoning, we have constructed 
a statistical model that captures the essential physics and can be solved exactly. The model 
has two principal ingredients. First, the unperturbed hamiltonian H is replaced by an 
ensemble of hamiltonians, each member of which has the same eigenenergies and eigenstates 
as Hq, but which eigenstate corresponds to which eigenvalue is chosen at random. It is 
convenient to think of these random assignments as permuting the eigenvalues, while the 
eigenstates remain fixed; so each permutation a, which sends the iV-tuple (1, 2, . . . , N) to 
(er(l), cr(2), . . . , a(N)), says that the jth largest eigenvalue corresponds to the eigenstate 
that in the original hamiltonian corresponded to the <7 _1 (j)th largest eigenvalue. Second, 
the perturbation V is replaced by an ensemble of perturbations; the distribution of the 
perturbations (viewed as matrix elements between various eigenstates) is assumed to be 
invariant under the permutations a of the eigenstates. We should note that we are permuting 
only the unperturbed eigenstates whose unperturbed eigenvalues lie within a microcanonical 
energy window Wmc(E, AE) centered at the mean energy E of the initial state. Also, the 
perturbations V are truncated so that they couple only the eigenstates within that energy 
window. The following result then follows (see the Supplementary Methods for the details 
of the derivation): 

7«\ ^ (N- N pc ) A° +N (N pc - 1) A MC 
h,v N pc (N-l) 

(2) 

Here N gives the number of the states in the microcanonical window Wmc{E, AE), and 

(...) stands for an average over the uniformly distributed permutations a and the 
\ / er,V 

permutation-invariant-distributed perturbations V. It will be convenient to introduce, for 
any quantity Ba that depends on the values of the quantum numbers n, the microcanonical 
average as 

MC E [B H ] = 1 J2 B * 

n:E a £W MC (E,AE) 
8 



where N is the number of unperturbed eigenstates whose eigenvalues are in the microcanon- 
ical energy window. Then the typical number of the (interacting) principal components is 
given by iV pc = l/^t? c{ . } 0, V)j . where Vuc{n}( a ^ ^) - MG e vfH^, V) is the micro- 
canonical average of the inverse participation ratio of the noninteracting eigenstates over the 
interacting ones (in other words, 7/1 = ^2 a |(a|n)| 4 ], A = (^mit.|^4|^init.) is the value of A 
in the initial state (it is written as an expectation value even though, given our choices, the 
initial state is also an eigenstate of A), and A MC = MCe [A^] is the microcanonical average 
for the observable. 

What we suggest now is to treat the "real-world" hamiltonian if as a single instance, 
(a = cr rcal - world 5 v = \/ real - world ) 5 f the ensemble of hamiltonians considered above. In this 
case, the ensemble average on the left-hand side of equation (J2J) will serve as the "best 

oo 

predictor" for the "real-world" value of the infinite-time average of the observable, A ■ 
Notice that, conversely, one of the constituents on the right-hand side of equation (|2j) will, 



in its turn, also have to be estimated by a "best predictor." Namely, while A and Amc are 
the same for all (a, V r )-instances and thus can be extracted from a single instance (which 
can be the "real- world" instance), the quantity Vuctnyi^^ ^0 ( on which N pc depends) varies 
from instance to instance and so its averaging over instances has a nontrivial effect. Since, 
however, we only have access to a single realization of H , we must use the "real-world" value 

t *■ J 1/t- l/l oo "i~rio hocf ocfiTYi o4-o Tr\v 4" ri a on ciqvyi r\ I o fnroro rro- / nr\ *■ ' 



of VucM\( a -> ^) as ^he "best estimate" for the ensemble average: \VMcm\( a 'i V) 



is the "real-world" 



„W l -.real- world f7rcal-world\ ^.i®} .,,L„ rQ „{&} A /IT 1 r,{ a } 

^MC{n}\ a > V > ~ VMC{n}i Wnere V MC{n} ~ ML E \Vft 

microcanonical average of VMCuyi 17 ^ ^)- 

An undesirable feature of the expression in equation (j2j) is that it depends on the number 
of states in the microcanonical window, N. However, in the limit where N greatly exceeds 
the number of principal components, iV ^> N pc , N disappears from the expression. We 
finally obtain the desired relationship between the infinite-time average and the initial value 
of our observable of interest: 

A°°= VMC{n} A° + (l - ^ { MC{n}) A ^ (3) 

This is almost our original estimate based on intuitive arguments, except that instead of 
the typical value of the IPR of the perturbed over the unperturbed states, r)a , we have 
the microcanonical average of the IPR of the unperturbed states over the perturbed states. 
Note the high degree of universality: the same parameter Vucm * s use d regardless of what 
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FI(jr. 2. Testing the predictive power of equation ([3)l. Our integrable system is a single particle on a two-dimensional 
lattice with an Aharonov-Bohm flux. Lattice parameters are the same as for Figure [T] Each row of plots corresponds to a 
different type of non-integrable perturbation, as follows. a,d, A real Gaussian random matrix perturbation acting between the 
eigenstates of the unperturbed lattice. b,e, A single impurity of fixed strength. c,f, An Anderson-type disorder. In plots a,d, 
and c,f, the data points marked with (red) squares are for the case e = 0.25; with (blue) circles, e = 0.5; and with (green) 
triangles, e = 1. In plots b,e, the (red) squares are also for e = 0.25, while the (green) triangles are for e = 2 X 10 5 . Plots 
a,b,c show the infinite time average, obtained from exact time dynamics, of the "equipartition measure" E x — E y versus the 
value this measure had in the initial state. The straight solid lines are the predictions of equation J3J. Included as insets are 
the representations of the lattice and its perturbation. In plots d,e,f we show the level-spacing histograms for the perturbed 
hamiltonians, which show how close the systems presented in plots a,b,c are to integrability or to well-developed chaos. The 
horizontal axis, Ae, is the level spacing in the unfolded spectrum. The closer the actual distribution to the curve representing 
the Poissonian level statistics, the more integrable the system; the closer the actual distribution to the curve representing the 
Gaussian Orthogonal Ensemble (GOE) or Seba or the Gaussian Unitary Ensemble (GUE), the more chaotic the system. See 
the main text for details. 



observable A one is interested in (which means — since A is diagonal in the unperturbed 
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basis — regardless of which integral of motion of H (n) one is interested in), or with which 
initial state one starts from among the unperturbed eigenstates. 

Testing the analytical prediction The relationships in equations (T5]) and ([3]) are our cen- 
tral results; we now test them against exact time dynamics of particular physical systems. 
Our integrable system will always be a 33 x 33- site two-dimensional lattice. We also assume 



that 
each 



)oth x- and ^/-cycles of the lattice are threaded by an Aharonov-Bohm (A-B) solenoid 



31 



33j . If the lattice is imagined to cover the surface of a torus, the flux is produced by 
two solenoids: one toroidal, contained within the torus of the lattice, and one straight, pass- 
ing through the hole of the torus. The corresponding magnetic field fluxes are assumed to be 
highly irrational but weak (i.e. of the order of one) multiples of the elementary magnetic flux 
quantum; they are also presumed to be mutually irrational — see the first subsection of the 
Methods section, below. Their purpose is to destroy the dihedral symmetries of the square 
lattice while preserving the conservation of momentum in both x- and ^/-directions. This lifts 
the degeneracies in the unperturbed spectrum and randomizes the sequence of appearance 
of the quantum numbers (two components of the momentum vector) along the energy axis. 
We investigated three different types of the integrability-breaking perturbation. The first 
perturbation (see Figure Eh) completes the lattice to a single instance of a Deformed (real) 



Gaussian Random Matrix Model 19j] that couples the eigenstates of the unperturbed lattice. 



The second perturbation (Figure [2b Ws a single impurity of a fixed strength. This system 



is a lattice version of Seba billiards |34j] „ which are known to lie in between integrable and 
completely quantum-chaotic systems [34 L The third example (Figure [2b) is a single instance 
of the conventional Anderson disorder 211 ] . with a rectangular distribution of the strength 
in each of the impurities. (We should mention that ergodicity in the context of an Anderson 



lattice was also studied in ref. (35|.) While the first and the second examples correspond to 
permutation- invariant perturbations, the third one does not. In particular, in comparison to 
the first two types of perturbation, in the Anderson case there are much stronger correlations 
between certain kinds of matrix elements. Namely, for each matrix element, consider the 
momentum difference between the states that the matrix element connects. If one looks at 
the set of matrix elements which have the same momentum difference, one finds that from 
one realization of the perturbation to the next, their values all change by the same factor. 
The observable of interest is the "equipartition measure" — the difference between the x- and 
y-hopping energies, E x and E y — whose value in a state of a thermal equilibrium is always 

11 



zero, thanks to the x -H- y symmetry. The infinite-time average of the "equipartition mea- 
sure" E x — E y is shown as a function of its initial value. As the governing parameter e we 
have chosen the ratio between the root mean square of the off-diagonal matrix elements V^v 
(the same for all the pairs (n, n')) and the typical energy spacing at the energy of interest: 



e = Vqp(E), where V — yl^nn'l 2 and p(E) is the density of states at the energy E. 

For each type of perturbation, we time-evolved 201 different initial states, taken to be 
all the eigenstates of the unperturbed lattice whose eigenvalues came from a representative 
microcanonical energy window; the middle eigenstate (the 101st one in the order of increasing 
energy within the window) had the energy of —1.5 J, where J is the hopping constant 
(see the first subsection of the Methods section). A single instance of the corresponding 
hamiltonian (which was random for the cases a and c) was used in all three cases. The 



different initial states have different values for (E x — E y ) , lying between some minimum 

n 

and maximum values; we needed to group the initial states into sets with similar (E x — E, 



yj 

-0 



values. Thus we partitioned the interval from the minimum to the maximum (E x — E y ) 

o 

value into subintervals centered at —3 J, —2.5 J, . . . 3 J, each of width A (E x — E y ) = .5 J. 



All the initial states whose (E x — E y ) -values fell into the same subinterval constituted a 



"group of initial states with similar (E x — E y ) -values." The points shown in plots a,b,c 
correspond to the groups: the x-value of a point is the center of the subinterval defining 

1 oc 

the group, while the y-value is the group average of (E x — E y ) . The theoretical curves, 



produced by equation (jHJ), also involved an averaging: the value of (E x — E y ) which is fed 



into equation (J3J) was, for each group, the group average of (E x — E y ) rather than the center 
of the subinterval; this explains why the theoretical curves are not exactly straight lines. 
As Figures |2^,b,c show, the prediction of equation fl3]) agrees very well with the numerical 
results. 

Figures |2H,e,f are there to show where on the continuum between integrability and well- 
developed chaos our various systems lie. The discrete points are the plots of the actual level 



spacing statistics of our systems, properly unfolded |36| ; for the cases d and f, we also average 
over 16 realizations of the perturbation V. For comparison, we also plot the curves corre- 
sponding to the completely integrable and (the appropriate) completely quantum-chaotic 
systems. The curves labeled "Poisson" arise for Poisson level statistics, which is charac- 
teristic for integrable systems. The most recognizable feature is the nonzero value at zero 
spacing, representing the absence of level repulsion in integrable systems. The Gaussian 
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Orthogonal Ensemble (GOE) curve is valid for systems with well- developed quantum chaos 



in the presence of time-reversal invariance; the Seba distribution 34J holds for systems with 
singular perturbations; and the Gaussian Unitary Ensemble (GUE) is used in the cases of 
quantum chaos without time-reversal invariance [36j. We see that when e = 0.25, the level- 
spacing statistics of our systems is intermediate between those of the integrable-like and the 
appropriate completely quantum-chaotic-like distributions. As e increases, the level spacing 
distributions for the real Gaussian, singular, and Anderson perturbations converge, respec- 
tively, to the GOE, Seba, and GUE predictions. (The reason for the GUE statistics in the 



Anderson case is the Aharonov-Bohm flux, which breaks the time-reversal invariance 32 ] . 
In the case of the first model the statistics remains of a GOE type, since the perturbation 
matrix elements were artificially fixed to real values.) The expression in equation Q is thus 
confirmed in the full range from the integrable regime all the way to the well-developed 
quantum chaos. 

To relate our predictions to the system parameters, we further connect the inverse par- 
ticipation ratio Vyicm — otherwise emprirically irrelevant — to the governing parameter e. 
This allows us to trace the integrability-to-chaos transition, i.e. express the memory of the 
initial state through the strength of the non-integrable perturbation. Figure [3] shows, for a 



particular group of the initial states, the values of the "equipartition measure" (E x — E y ) 
as a function of the parameter e. The numerical values are drawn from the set used to 
produce Figure [2j To obtain the theoretical prediction, we tabulate numerically the inverse 
participation ratio r/i" , averaged over the unperturbed states |ao), for a Deformed (com- 
plex) Gaussian Random Matrix Model. For large values of e, we also use known theoretical 
results |37|, |38[ for the so-called strength function, |(a |a)| , combined with the assumption 



391 ] of the Gaussian character of the fluctuation of (ao\a): 



— oo n — o 

A -A MC « T-y^U -A MC ) , (4) 

where q = 3 for a non-integrable perturbation that belongs to the Gaussian Orthogonal 
class, and q = 2 in the Gaussian Unitary case. Note that the Aharonov-Bohm flux present 
in the example of the Figure [3] implies the latter. Note also that the result in equation (j4j) 



for the Gaussian Orthogonal case can be confirmed directly 38(, without the assumption of 
Gaussianity of (a \a). The strength function for an equidistant spectrum perturbed by a 
real matrix whose matrix elements have the same magnitude and random signs — the case 
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closely related to a Poisson spectrum perturbed by a real Gaussian matrix — first appears 
in the classic Wiener's paper J40(. The strength function for the latter system per se is 



predicted in ref. 



411 ]. (See the second subsection of the Methods section for more details.) 




FIG. 3. The memory of the initial state as a function of the strength of the 

integrability-breaking perturbation. For the case of Anderson disorder, we compute, for three 
similar initial states, IV'init.) = \n x = +7, n y = 0,±l), the infinite time average of the "equipartition 
measure" E x — E y as a function of the disorder strength e. The red crosses correspond to the 
averages over the this group of states. The red solid line corresponds to the formula (|3|) combined 



with a numerically tabulated values of the inverse participation ratio r\ 



M 

MC{n} 



for the case of the 



Deformed (complex) Gaussian Random Matrix Model. Finally, the black dotted line is given by 
the asymptotic analytic prediction in equation (JH). The lattice parameters are the same as for 
Figure [TJ 

DISCUSSION 



We have thus demonstrated that, in the case when the integrability-breaking perturbation 
obeys no selection rules and when the spectrum of the underlying integrable system is 
"sufficiently irrational," it is possible to characterize the memory of the initial values of 
the unperturbed integrals of motion, as the strength of the perturbation is increased, by 
a simple and universal expression. The expression was verified for three different types 
of perturbations away from integrability, including a Deformed Gaussian Random Matrix 
Model, an isolated impurity, and the two-dimensional Anderson model; in each case, the 
expression works for a full range of perturbation strengths from completely integrable to 
completely chaotic. 
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Our predictions should be testable experimentally with techniques that are already avail- 
able, or nearly so; likely contexts include those of the investigations of Anderson localization 
with cold gases in optical lattices 22|, |23J], or those of quantum mirage configurations of a 
quantum corral [42], modified to a Seba-type billiard. In general, memory effects should be 
visible as soon as the discreteness of energy levels is discernible. 

Future work must try to come to terms with cases where the partial preservation of the 
integrals of motion is enhanced via some selection rules obeyed by the perturbing potential. 
The two most empirically relevant examples of selection rules come from the limitations 
imposed by the few-site nature of the hopping terms in typical lattice hamiltonians, and 
by the few-body nature of interactions in many-body systems. Encouragingly, some work 



relevant to these cases has already been done; ref. 43| studied the role of the former type 
of selection rules, while the effects of the latter type on the structure of the eigenstates 
has been investigated in a number of works; for a review, see ref. [19| and the references 



therein, particularly ref. 44j . More generally, one should start systematically increasing the 



complexity of the topology of the network of transitions. 



METHODS 

Numerical models used to verify equation |3|). As an example of an unperturbed hamil- 
tonian H , we use a N x x N y (33 x 33 in the numerical examples considered) two-dimensional 
lattice with periodic boundary conditions and odd N x , N y : 

N X -1 | Ny-1 

2 2 

Ho = ~J E E (e~^l(;* + l,iy))<(i*,i,)l + 

j^- N X -1 Ny-1 

. 27vA.riy 



e »v \(j x ,j y + l))((j x ,j y )\+h.c.J, (5) 

where \(j x , j y )) are the eigenstates of position. Both x- and y-hopping constants have the 
same amplitude J. We also assume that both x- and ?/-cycles of the lattice are threaded by an 
Aharonov-Bohm (A-B) solenoid each. This leads to the hopping constants acquiring complex 
phase factors, with phases —2irAn x /N x and —27rAn y /N y respectively. Here An x = 4> x /<j>o 
(An y = 4> y /4>o) where <p x ((f> y ) is the magnetic flux through the x-cycle (y-cycle), 0o — 2nhc/q 
is the elementary quantum of magnetic flux, c is the speed of light, and q is the electric charge 
of the lattice particle. The flux is assumed to be very weak, An x (An y ) ~ 1. Note also that 
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after a suitable gauge transformation, the complex hoppings can be replaced by real ones, 
supplemented by twisted boundary conditions J3lM33| . The eigenstates of the unperturbed 
hamiltonian are plane waves \(n x , n y )), 

((j x , j y )\(n x , n y )) = (N x N y r^e Xp {^^ + .^ 

The linear momentum quantum numbers n x and n y (\n x \ < (N x —l)/2 and \n y \ < (N y — l)/2) 
constitute a set of the integrals of motion. Eigenenergies of the unperturbed hamiltonian 
are 

E nx , ny = -2J {cos[2-K(n x + An x )/N x ] + cos[2-K(n y + An y )/N y }} . 

The purpose of introducing the A-B flux is to generate an equi-energy surface E nxjUy = E 
that is sufficiently irrational with respect to the lattice of the integer quantum numbers n x 
and n y . To this end, the "defects" An x and An y must be both irrational and mutually 
irrational (see the first subsection of the first section of the Supplementary Methods and 
Supplementary Figure S2). In our numerical calculations, we used An x = ip/4 and An y = 
e/10, where if = 1.618 ... is the golden ratio and e = 2.718 ... is the base of the natural 
logarithm. 

As the first example of an integrability-breaking perturbation V, we considered a member 
of a Gaussian Orthogonal Ensemble acting between the eigenstates of the unperturbed lattice 
(Figure^): 



Vfif{' — V 1 + (WVo£nn' ' 

where the ^n' are N x N y (N x N y + l)/2 independent real Gaussian-distributed random vari- 
ables of unit variance and zero mean, with the remaining N x N y (N x N y — 1)/2 matrix elements 
controlling the hermiticity of V. 

The second example (Figure [2b) is a singular perturbation, 

V(J X , j y ) = V N x N y 5 {jxJy)A o,o) , 
which produces a matrix with all-equal matrix elements: 

Vhh' = Vo ■ 
In this case, the eigenstates of the perturbed hamiltonian can be found exactly. They read 

| N X -1 | Ny-1 

IK, n y )) 



i«> = cr s £ E 



N x -i Ny-l a n x ,n y 

2 n y— 2 — 
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where the corresponding eigenenergies E a are the solutions of the algebraic equation 



1 N X -1 . Ny-1 






_ if 1 JV i « 
n x - — n y = -?— 


1 

- E 

iT>X i Toy 


1 

~V 



and the normalization factor is defined by 

| N x -1 | ^-1 

(cro~ 2 = E E (^.--Ewj" 



A solution of this form was at first obtained by Seba 20J], for a flat continuous billiard. A 
similar problem involving a two-dimensional lattice with periodic boundary conditions in 
one direction and a trapping potential in another was recently analyzed [45|. 
Finally, we consider an Anderson-type disorder (Figure [2b): 

where the Cj x ,j v are N x N y real independent variables, distributed uniformly between —1 and 
+1; W is the Anderson disorder parameter. The interaction strength parameter Vq used in 
the previous cases corresponds to the r.m.s. of the (generally complex) off-diagonal matrix 



elements V nn >, V = ^\V m ,\ 2 = W/(2^^NjQ. 

Relevant results on the IPR in random matrix models. To obtain the solid line in Fig- 
ure El we computed the inverse participation ratio r/a" , averaged over all unperturbed states 
|o!o)) f° r a N x N Deformed Gaussian Random Matrix h = h + v, with iV = 2000. The 
"integrable" part of the matrix, ho, was represented by a diagonal matrix whose iV diagonal 
entries were given by N independent real random numbers uniformly distributed in the in- 
terval [— (N/p)/2, +(N/p)/2]. Here, p is the density of states. The "non-integrable part", v, 
was a random matrix drawn from the Gaussian Unitary Ensemble of random Hermitian ma- 
trices: (real) diagonal matrix elements and the real and imaginary parts of the off-diagonal 
matrix elements above the diagonal (those below the diagonal then being fixed by her- 
miticity) were given by independent, Gaussian-distributed random numbers with standard 
deviations o"di ag . = Vq and cr fj-_di a g., Ro = Coff-diag., i m = Vo/\/2, respectively. For completeness, 
the Gaussian Orthogonal case would give <7di ag . = \/2Vo and cr ff-dia g . = Vq. In both cases, Vq 



fixes the mean square of the off-diagonal matrix elements: \v a a ' | 2 = V 2 . 
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Equation (J5J) uses the following asymptotic expression for the inverse participation ratio 
//: 

(6) 



' e>i 2vr% 2 ' 
where q = 3 for a non-integrable perturbation that belongs to the Gaussian Orthogonal 



class, and q = 2 in the Gaussian Unitary case. In refs. 37], |38|, |40) and 41] it was shown 
that for e ^> 1, the strength function converges to 

i nr/2 

e >i ~p {E ao - E a y + (hT/2)< 



K«o|«)| 2 « —tb F ^, 7^77^ » (7) 



in both Gaussian Orthogonal and Gaussian Unitary cases. Here, V is determined by the 



Fermi Golden Rule: V = l-nV^plh. Also, it is known 39[ that in the same limit, the 
individual coefficients (ao|«) behave as independent Gaussian random variables, with zero 
mean and with a standard deviation governed by the expression (JTj). Recall that the fourth 



moment of a Gaussian distribution of a variable £ is related to the second one as 46] 
l£| 4 = Q ( l£| 2 ) • Also, it is important that when e ^> 1, the energy width V in equation (EJ) 
contains many levels, and thus the sum in the definition of the inverse participation ratio, 
r)s% = ^2 a |(q;|q;o)| 4 , can be replaced by an integral. The asymptotic formula (JHJ) for the 
inverse participation ratio immediately follows. 

Ab initio computations. All numerical time evolution was computed from full exact 
diagonalization of the hamiltonians using the Mathematica computer package. 
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SUPPLEMENTARY FIGURES 
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Supplementary Figure S4. The energy-ordered sequence of the eigenstates of a "real- 
world" integrable system viewed as a random permutation of an analogous sequence 
of the eigenstates of a fictitious "reference" hamiltonian. 
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Supplementary Figure S5. A test verifying that the integrable hamiltonian in equa- 
tion (5) of the main text has a "sufficiently irrational" equi-energy surface. We find 
no correlations between the quantum numbers of pairs of eigenstates that are nearest neighbors 
in a sequence ordered by increasing energy: a, A two-dimensional histogram for the distribution 
of the ((n x ) ao , (n x ) ao+ i) pairs of the values of the x-momentum of two consecutive eigenstates 
(the index «o labels the eigenstates in order of increasing energy), showing no apparent correlation 
between (n x ) aQ +i and (n x ) ao . b,c, The same histogram as in plot a for two fictitious sets of eigen- 
states obtained via two random permutations of the original sequence (see also Supplementary 
Figure IS4|) . presented for comparison with the histogram in plot a: the latter, which is the one for 
the actual system, looks as random as the histograms in plots b and c, which were generated by 
random permutations of eigenstates. The integrable hamiltonian was a 33 x 33-site lattice pierced 
by a weak Aharonov-Bohm (A-B) flux. The histograms include states from a representative micro- 
canonical energy window comprising 201 eigenstates with the energy of the middle (101st) state 
being E = —1.5 J, where J is the hopping constant. Individual boxes of the histogram contain 
3x3 = 9 values of ((n x ) ao , (na;) Qo +i). For other details, see the beginning of the first subsection 
of the Methods section of the main text. 
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SUPPLEMENTARY METHODS 

A model supporting the formula in equation (2) of the main text 

The unperturbed integrable hamiltonian. 

Consider the hamiltonian of an integrable quantum system under a nonintegrable per- 
turbation V: 

H = H (k) + V , 

where n = {fii, fi2, ■ ■ ■ , fid} is a complete set of the integrals of motion of the hamiltonian 
Ho(n), and d is the number of the degrees of freedom. If the system is large enough, one 
should be able to identify a macroscopic energy scale S that governs the regular (as opposed 
to fluctuating) part of the dependence of thermal expectation values of observables on energy. 
Under these circumstances, it is possible to choose, for a given energy E, a microcanonical 
window 

W MC (£, AE) = [E- AE/2, E + AE/2] (S-8) 

that is small compared to the macroscopic scale £ but still large enough to contain a large 
number of the eigenstates. Let us arrange the eigenstates of Ho(n) inside the window in the 
order of increasing energy and label them using an index a,® — 1, . . . , N: 

a > a' => E n > E n (S-9) 

a Q 

E nao=1 *E-AE/2; E Rao=N » E + AE/2 , 

where N is the number of the eigenstates the window contains. The hamiltonian of the 
unperturbed system reads 

N 

H =^2 E cco\n ao ){n ao \ + • • • , (S-10) 

Q0 = 1 

where \n) and E^ are the eigenstates and eigenenergies respectively, 

E — E- 

and ". . ." stands for the states outside the microcanonical window. 
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Now, we suggest replacing the hamiltonian in equation ( IS-lOj) by an ensemble of random 
hamiltonians: 

N 

Ho(&) = JZ E<x \n ao (cr)){n ao (cr)\ + ... (S-ll) 

«o=l 

where n„ is a "reference" sequence of the eigenstates, and the permutations a are assumed 
to occur with equal probability p = 1/N\ . The "real- world" sequence n ao is supposed to 
be represented by a particular realization of a random hamiltonian from equation (IS-llj) 
corresponding to a particular permutation <j real - world : 

Notice that the spectra and the sets of the eigenstates are identical for all members of the 
ensemble; however, the order in which the eigenstates appear on the energy axis is different, 
and is realized randomly The relationship between the reference sequence of eigenstates 
and the sequence of eigenstates as it is for the actual, "real-world" integrable hamiltonian 
is illustrated in Supplementary Figure [SU 

The replacement of the actual integrable hamiltonian by an ensemble of random integrable 
hamiltonians is justified if the sequence of eigenstates in the actual integrable hamiltonian 
is indeed "sufficiently random." In the main text, this property was phrased in terms of the 
integrable hamiltonian having a "sufficiently irrational" equi-energy surface: let {{nk)v) be 
the sequence (labeled by v) of the values of the /cth quantum number as one is going from 
one eigenstate of H (n) to the next in the order of increasing energy. The equi-energy surface 
is said to be sufficiently irrational if, for every k, the sequence {{fik)v} passes any simple 
statistical test for randomness. One such test is explained in Supplementary Figure [S5] on 
the example of the integrable hamiltonian defined in detail in the first subsection of the 
Methods section of the main text. 

The non-integrable perturbation. 

Let us first assume that the perturbation V is weak as compared to the width of the 
microcanonical window, 
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so that it mostly couples the states inside the microcanonical window of equation ( IS-8J) to 



themselves. In this case we can approximate the perturbation by its truncated version, i.e. we 
neglect all the matrix elements that involve the states outside the window in equation (1S-8[) : 

V m > ^ =}► En e W MC {E, AE) & E ni e W MC {E, AE) , 

Let's further assume that the perturbation V does not obey any apparent selection rules. 
In this case, V can be thought as a particular realization of a random matrix, with the dis- 
tribution of the matrix elements, w(V), invariant under the permutations of the eigenstates: 

w(U a VU; 1 )=w(V) , (S-12) 

where 

The Gaussian Orthogonal (respectively, Unitary) ensemble is the most obvious candidate for 
such a distribution: its distribution of the matrix elements (real for the Gaussian Orthogonal 
ensemble and complex for the Gaussian Unitary), 

^goe(gue)(U) cc exp[-Tr(U 2 )/4U 2 ] , 

is invariant under orthogonal (respectively, unitary) transformations, and in particular under 



permutations. Here, Vq = y \VftH'\ 2 is the standard deviation of the off-diagonal elements. 
A (fixed, non-fluctuating) singular perturbation, 

Wsingular(VO °C J| |^[ 5(Vftfi> ~ V ) , 

n n' 

is another example of a "permutation-invariant" distribution of the perturbation matrices. 

A proof of the formula in equation (2) in the main text 

For a given permutation a and a given set of the matrix elements of the perturbation V, 
the infinite time average of the quantum expectation value of an observable A will be given 
by 

A°° (a, V) = lim - [ dt(ip(t)\A\ij){t)) 

T^OO T J Q 

= ^2\(a\n init )\ 2 (a\A\a) 



7 CJ (T '^ + E Fftft in]t (cr,V)Aft, (S-13) 






2(3 



The significance of singling-out the first term in the last line will become apparent later in 
the derivation (see the remark after equation (IS-221) ). Here 



A = A H . nit 

is the initial quantum expectation value of the observable; \a) is an eigenstate of the per- 
turbed hamiltonian, of an eigenenergy E a : 

H\a) = E a \a) ; 

the quantity 



a=l 

is the inverse participation ratio of an unperturbed eigenstate \n) over the perturbed eigen- 
states | a); the function 

N 

Fft } (tT,y)^E|(a|n)| 2 |Hfi'}| 2 (S-15) 
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is the so-called F-function . Here and below, we assume that the initial state, \ip(t — 
0)) = l^init.), is one of the eigenstates of the unperturbed hamiltonian from equation (IS-lip 
that belong to the window of equation (1S-8J) . and that the observable A is diagonal in the 
eigenbasis of the unperturbed hamiltonian: (n\A\n') = A^Snn 1 - As it was mentioned above, 
we neglect all matrix elements of V that couple the inside of the window with the outside. 
One can replace the n indices in equations (1S-14HS-15]) by the a indices: 

»£V, v) s „£> w ( ff , v) (s-i6) 

*&(*. V) = f*'^ w («r, V) , (S-17) 

see equation (1S-9|) . The infinite-time average of the observable now becomes 

where ao.initX* 7 ) — ^(^o.init.) is the position of the initial state in the sequence given by a 
permutation cr, and a:o,mit. is the position of the initial state in the "reference" sequence: 



n «0,init. ~~ ninit - 
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Next, consider the average of the infinite-time average of the observable over the realiza- 
tions of the perturbation V: 

Q 0^a ,init.( CT ) 

Notice that we used the property that, by construction, the IPR and the F-function, if 
expressed as functions of a , depend on the permutation a only through the perturbation 
matrix elements V a a ' = V^ a ( CT )^ , (a)- Therefore, by virtue of the invariance of the distribu- 
tion of V with respect to the permutations, (see (1S-12J) ). the average IPR and the F- function 
(again, expressed through the a indices) do not depend on the permutation. That is, for 
any two permutations a and a', 

F {a \) (<j) = (F {a \) (a') = (F {a \) . (S-19) 

Now, let us consider the average ( A ) , which is the average of the infinite-time 

\ IV, a 

average of the observable over both the ensemble of V°s and the ensemble of the unperturbed 



hamiltonians whose members are parametrized, as before, by the permutations a. Since A is 
invariant over permutations a, the cr-average of the first term in the r.h.s. of equation (1S-18J) 
becomes 

» \ T"\ _/„(«> ' — 



Ch^ = W^).^ • < s - 2 °) 



where 



0,init.W/V la \'™^t»/ v ' ' / (T,V 



vVc { ft } (°, V) = iV- 1 J2 <\°> V) (S-21) 



00=1 

N 



ft:E a eW MC (E,AE) 

is the microcanonical average of the inverse participation ratio, for a hamiltonian parametrized 
by a permutation a and the perturbation matrix V. Note that formally speaking, the o~- 

averaging in the r.h.s. of equation (JS-20P is not necessary. However, in any approximate 

>} 



V^-average over a finite series, the a-dependence of \Vucfn,}( a ' V) 7 - wn ^ inevitably appear 
again. 
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The cr-averaging of the second term requires a two-stage procedure. First, we will average 
only over the permutations for which the initial state index ao.init. (0") is fixed to a particular 
value ao.init.- Using the property flS-191) we get 



o.imt.v 1 1 v ° I <t|<* init »=a* 

a o^a ,init.(°") 



T ( F{a} n x ^ 

E (*%_,.)» xi^-^-D-'U'-iMo)} (s-22) 



u ~ O.mit. 



where 



iV 



A MC = N- 1 J2 A n 

H:Efi£W MC (E,AE) 

is the microcanonical expectation value of the observable. Notice that Amc depends neither 
on permutation a nor on perturbation V. Note that if we were to retain the a® = a$ init 
term on the left hand side of the equation ( 1S-22I) . the factorization into a product of an 
intitial-state-dependent and intitial-state-independent factors on the right hand side would 
not be possible; this, in turn, would make the subsequent simplifications impossible as well. 
This justifies the singling out the first term in Eqn. C IS-13J) . which is the term responsible for 
the memory of the initial conditions. 

The second factor in the r.h.s. of equation ( 1S-22J) does not depend on the initial state 



index a$ init at all. Averaging of the first factor over a^ init gives 



£ «"4 ),).,. =i-(^W(<^)K ■ (s-23) 



■o"o,init./v7a* .. \ ■ ™"V*J < ' / CT , y 



o.mit. 



Combining the particular results in equations ( 1S-20IIS-22IIS-23I) . we get 



M 



A )-=^o^y)),A + 



-0 



1 - (v { $ C{ n } & V)) J x [A MC - (AT - 1)-\ A ° -A MC )} (S-24) 



,v / 

(N - N pc ) A° +N(N pc - 1) Amc 
A pc (A-l) 

where the number of the principal components is given by 

1 



iV P c - {«} 



V MC{n} 

This concludes the proof for the formula in equation (2) of the main text. 
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